%UPA导向矢量函数
function y = UPA_fun(phi,theta, Nx, Ny)
for m= 0:Nx-1
    for n= 0:Ny-1
        y(m*(Ny)+n+1) = exp( 1i* pi* ( m*sin(phi)*sin(theta) + n*cos(theta) ) );
    end
end
y = y.'/sqrt(Nx * Ny);
end
